Analytic theory of optical nano-plasmonic metamaterials 
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Recent advances in nano-fabrication techniques allow for the manufacture of optical metamateri- 
als, bringing their unique and extra-ordinary properties to the visible regime and beyond. However, 
an analytical description of optical nano-plasmonic metamaterials is challenging due to the char- 
acteristic optical behaviour of metals. Here we present an analytical theory that allows to bring 
established microwave metamaterials models to optical wavelengths. This method is implemented 
for nano-scaled plasmonic wire-mesh and tri-helical metamaterials, and we obtain an accurate pre- 
diction for their dispersive behaviour at optical and near-IR wavelengths. 



The metamaterial paradigm 'function from structure' 
has in recent years opened the door to exciting new mate- 
rials with unprecedented properties such as negative re- 
fractive index Moreover, in combination with trans- 
formation optics, it allowed a conception of dramatic 
new functionalities such as perfect lensing cloak- 
ing [H 0| or broadband slow and stopped light Q ■ As 
the meta-atoms that make up a metamaterial charac- 
teristically consist of metal structures that are smaller 
than the wavelength, metamaterial principles such as in- 
visibility cloaking 0, @| were first experimentally demon- 
strated in the microwave regime. Recent advances in 
nano-fabrication methods, particularly bottom-up and 
self-assembly techniques now allow for the design and 
manufacture of artificial nano-scaled metallic structures 
with complex structural features [H-Q forming nano- 
scaled metamaterials operating at optical (near infrared 
and visible) frequencies. 

In the conception of metamaterials, analytical theories 
have been the trailblazers for their design and opened 
our understanding of their physical properties. For chi- 
ral metamaterials this has been particularly successful 
in the microwave region At optical wavelengths, 

the forces excerted by incident lightficlds on the nearly- 
free (conduction) electrons of a metal lead to collective 
electron-photon oscillations that typically occur on scales 
of just a few tens of nanometres or less. These collec- 
tive electron-photon oscillations can either be localised 
(localised surface plasmons, LSPs) or propagating (sur- 
face plasmon polaritons, SPPs), forming areas of high 
intensity around and inside the plasmonic nanoparticles, 
the meta-atoms. In addition to nano-localisation of inci- 
dent light-waves, surface plasmons may efficiently allow 
for the reverse effect, i.e. the out-coupling to the far- 
field continuum of the near-field of emissive molecules 
and quantum dots placed adjacently to the metal meta- 
molecule, thereby acting as nanoantennas. Clearly, ana- 
lytical metamaterials models that are very appropriate 
for microwaves fall short of describing the optical re- 
sponse of nano-scaled metallic structures and nanoplas- 
monic metamaterials. Indeed, it was found through ex- 
perimental results and numerical calculations [HI [l2| 
that the actual performance of these structures is sig- 



Figure 1. A nano-wire made from metal of permittivity e m 
and surrounded by a host medium of en is completely pene- 
trated by an incident electric field E. The currents induced 
inside the wire give rise to a magnetic field \H\ shown at the 
cross-section surface of the wire. The amplitude of \H\ and 
the current density | jo | inside the wire are also displayed along 
a line parallel to the y-axis and passing though the centre of 
the wire (blue-shaded area represents the metallic wire). 



nificantly red-shifted compared to what one may have 
expected from their geometrical parameters 13j. 

In this Letter, we develop an analytic theory for the op- 
tical properties of metal nano- wires and propose a robust 
self-consistent method that can be used to bring well- 
established analytical metamaterials models from the mi- 
crowave domain to optical wavelengths. In particular, we 
will demonstrate the significance of our theory in brin ging 



the celebrated 'wire-mesh' metamaterial model [14j, [15 1 
to the optical frequency domain with a nano-scale wire- 
mesh and introduce an analytical model for tri-helical 
nanoplasmonic metamaterials. 

Let us start by focussing on a thin straight and 
infinitely extended metal nano-wire of radius r w 
and permittivity e m that is embedded in a dielec- 
tric host medium with permittivity Eh, as shown 
in figure [TJ The skin-depth S is dependent on 
the properties of the metal and given by: 6 = 

- ^(l - C))/((C) 2 + (I - e' m m 2/(kfr m ), 
where fcg is the wavenumbcr in free space and [i m the 
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Figure 2. The four solution (\zi\) plotted for Drude-silver 
metal, with the approximated solution (|z ap proa:|) of Novotny, 
which ignores metallic losses. 
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Figure 3. The H-function against the frequency for silver-Ag 
and gold-Au(o; mp = 1.3544el6(rad/s), g = 1.5367el4(l/s), 
£oo = 9.0685) in vacuum and eh = 2, and Aluminium-Al 
{io mp = 2.295el6(rad/s), g = 1.27el4(l/s), £oo = 0.8476). 



magnetic permeability of the metal. At microwave 
frequencies, any current induced on a metallic wire 
will flow on its surface, producing a magnetic effect 
around the wire and there would be no field inside the 
metal. Now, at optical regimes, the fields penetrate in 
the metal (since r w < S) and they induce an electron 
flow throughout the metallic nano-wire and not just on 
its surface. The current flowing at the centre of the 
nanowirc induces a magnetic field in a loop endorsing the 
current, as shown in figure [T] Of course, the magnitude 
of the induced magnetic field increases with the amount 
of metal it encloses and therefore it reaches a maximum 
value just outside the metallic nano-wire. Now these 
fields inside the metal slow down the effective/overall 
electron flow, which is modelled in this paper with the 
wavenumber 7 and significantly red-shift the effective 
response of the nano-wire. 

Hence,the wavevector describing the electron flow is 
kQy/s m and the effect of plasmon polarities (i.e. charge 
distribution inside the metal) is taken on board the 
wavevector 7, such that the overall field propa gati on in- 
side the wire is governed by Ki = yfcp£m — 7 2 [13|. The 
fields in the dielectric host medium are affected in a sim- 
ilar way by the plasmonic excitation, and therefore we 
have K2 = \J k^Eh — 7 2 - Hence, 7 dominates the red-shift 
to the dispersive behaviour of nano-metallic structures 
from the w frequency predicted from microwave analyt- 
ical models, to uj' that is actually observed from experi- 
mental results and numerical calculations. 

The relationship between the two frequencies (to and 
u/) is therefore dependent on 7 as [l3j |: 



in 



4r UJ qj / 
27TC 



an exact solution for (&0/7) and then we replace all ui 
parameters in established microwave analytical models 
for wire-mesh and tri- helical metamaterials with (JTJ), to 
obtain their optical dispersive behaviour. 

In order to find 7, we observe the field behaviour in 
figure [U which is identical to the TMq mode field profile 
for a cylindrical dielectric waveguide. The nano-wire thus 
effectively acts as a cylindrical dielectric waveguide since 
it allows the field to completely penetrate it. Hence, we 
use the waveguide theory to model the field penetration 
and plasmonic excitation for a metal nano-wire and we 
find 7 by solving : 



Kir w Jo(Kir w ) K 2 r w 



where J„ and are the cylindrical Bessel and Hankcl 
functions of the first kind. For thin wires (i.e. A 3> r w ) 
the Bessel and Hankel functions are expanded to first or- 
der in nr w terms and the above equation can be reduced 
to 0: 



(2) 



z 2 + a{z) + 



2e h (2a(z) - 2a{zf 
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(3) 



where a(z) = T — ln(2) + ln|z|, T = 0.577 is the Euler 
constant. We approximate ln|z| ~ 2^|^-^, which is 

valid for 0.5 < \z\ < 3. Therefore, a(z) ~ a + 2 (f^r)) 
where ao = T — ln(2). Also, where z (l3| | : 
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where the only unknown parameter is 7. Here, we find 



There are four solutions to (j3J) , plotted in figure [5] for 
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silver, which have a rather long and complex format: 
«o 
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and J = 3e m (2 + 5a + ag) - e h (5 + 6a + 2ag), G = 
2e m (a + 2)(a — 2) — eh(2aQ — 2a — 7) and H is a com- 
plex functions of Eh and e m , whose full format is shown 
in the supplementary material. In figurc[3J function H is 
plotted for different metals and host media. We observe 
that for frequencies where the real part of the metal's 
permittivity is (s m / £ h) < — 5> H is independent of the 
metal's and host medium's material properties. There- 
fore, to slightly simplify solutions in ([5]), H can safely 
be approximated by H = —1.537 + iO. 658712 for wave- 
lengths where (e m /£h) < —5. The four solutions of \z\ are 
plotted in figure [5] for a Drude-like silver of £oo = 4.028, 
plasma frequency uj mp = 1.39el6 (rad/s) and collision 
frequency g = 3.14el3(l/s). From figure [2] one can see 
that \z\\ and | ^3 1 are always degenerate. On the other 
hand, \z2\ and \z±\ are only degenerate for wavelengths 
where (e m /eh) > —5. Finally, we note that all our so- 
lutions are within the approximation limit 0.5 < \z\ < 3 
that we took earlier for ln\z\. 

Using equations ([T]) and (j4]), we can now obtain (fco/7) 
that is also plotted in figure [2] for all four z solutions. Fi- 
nally, the frequency shift occurring in dispersive metallic 
media, when they are completely penetrated by electro- 
magnetic waves is given by: 



(9) 



2ttc 



where u is the frequency found in any microwave analyti- 
cal dispersive model and a/ is the new frequency account- 
ing for the field penetration and plasmonic excitation of 
metals at optical frequencies. Note that z is dispersive 
(i.e. z{uj')) since it is a function of e m . Equation © 
can now be implemented easily into any established mi- 
crowave analytical metamaterial model, bringing them 
to optical wavelengths. Finally, ([9]) is valid for metal 
nano-structures of circular cross-section. If one wishes to 
obtain the relationship between u and u' for metal nano- 
structures of square or orthogonal cross-sections, then @ 
needs to be rewritten accordingly (and follow the same 



methodology). In this article, we choose to derive © 
for circular cross-sections, since at the nano-metre scale 
metal structures most commonly have approximately a 
circular-like cross-section. 

Novotny in [ijj], derived an approximate solution 
for (|5|) assuming loss-less Drude metals. His approximate 
solution Zapprox = ol\ + ol2~^j- , also plotted in figure [2] 
(blue dotted line) with our exact solutions. The approxi- 
mation provides an indication for the order of magnitude 
of the frequency shift for simple metallic designs, such as 
a dipolc optical antenna as it was intended to do (where 
L0p m is the plasma frequ ency of the metal and a\ and ui 



are coefficients given in [13j ) . 

We note that despite the fact that Novotny's approxi- 
mate solution [l3| gives quite reasonable results for sim- 
ple systems such as optical dipole antennas (rods of 
short length) , it cannot predict the properties of complex 
metallic geometries. In particular, metamaterials that 
quite often are complex designs of several wavelengths 
in size (i.e. wire-mesh, tri-hclical and gyroidal meta- 
materials). Furthermore, losses are an important aspect 
for metamaterial performance and cannot be neglected 



[17| . In contrary to 



our exact solutions in 



accounts for metallic losses that can significantly affect 
the performance of optical metamaterials. However, the 
most important advantage of our model is that it allows 
for any metal (not just Drude-metals) to be considered 
and in fact experimental data (in the form of a fitting 
function) to be used for e m to even account for inter sub- 
band transitions. Additionally, the model described in 
this paper can be also applied for both inhomogeneous 
and non-linear host media (i.e. gain and kerr media), 
by simply expressing Eh with the appropriate function. 
The derived analytical formulation in this paper unveils 
the physics governing the electromagnetic interaction of 
metal nano-structures with light in the infra-red and visi- 
ble regimes and consequently, modelling the homogenized 
electromagnetic behaviour of optical metamaterials. 

To demonstrate the generality of our theory on analytic 
methods, we now apply it to two different nanoplasmonic 
metamaterials, the wire-mesh Il4l Il5l Il8 | and the tri- 
helical metamaterial (THM) flpj . 

The classic wire-mesh metamaterial is composed of 
thin wires aligned with the three orthogonal axes (as 
shown in the inset o f fig ure IHb)) and behaves as an artifi- 
cial plasma 14j, [l_5|, [l8j . Its effective electric permittivity 
obeys the Drude model: 



1 



l"fUJ 



(10) 



where u p - mp is the cut-off frequency given by |l4 lla , ll8| : 

that 



co^V^l^y^T) + 0.5275)). Note 



p—wm 

we omitted eh in the above formula, since it is taken into 
account from our model in Now, if r w < S, we sub- 
stitute uj in (|10[) with §§§ . Assuming silver nano- wires of 
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Figure 4. (a)The analytical prediction (lines) and numerical 
calculations (points) for a silver wire-mesh with r w — 3.5nm 
and periodicity a — 40?im. (b) The cut-off frequency variation 
for different r w and Eh of perfect electric conductor (PEC) and 
silver wire-mesh nano-metamaterial. 



radius r w = 3.5ran arranged in a periodic square lattice 
of a = 40nm, we calculate e wm using (fTUj) and ([9]) and is 
plotted in figure Ufa) with numerical calculations [l9[ for 
comparison. The numerical values were retrieved from 
S-parameter calculations, using the method described 
in [2(J. We see a good agreement between the two over 
a large wavelength spectrum. It can be shown that the 
small discrepancies for higher Eh are due to the initial 
analytical model in (|10|) . and scale with the wavelength. 
Note that for z in ([9]), we used solution z±, since z\ 1 zi 
and Z3 gave unphysical electromagnetic properties for the 
optical wire-mesh metamaterial. In figure |U[b), we plot 
the analytical prediction (lines) of the wavelength where 
£wm = for various nano-wire radii (r w ) in different 
host media (eh) that show remarkable agreement with 
numerical calculations (points). It is noteworthy that 
the wire-mesh cut-off frequency of extremely thin wires 
significantly shifts from few tens of nanometre for PEC 
metal to few micrometres for silver. Also, small varia- 
tions on the optical wire radius produce extreme shifts 
for the cut-off frequency and our model predicts very ac- 
curately this dramatic behaviour. 

To further demonstrate the applicability of our method 
on any nanoplasmonic metamaterial, we now apply it 
to a chiral metamaterial. The tri-helical metamaterial 
(THM) was introduced in [l(| , as an analytic model for a 
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Figure 5. (a) The n_ and (b) n+ of a THM metamaterials 
of r w = 3.48nm, pitch p — 20nm, helix radius R — lOnm and 
periodicity a = 40nm. (c) The cut-off frequency variation for 
different r w and Eh of perfect electric conductor (PEC) and 
silver wire-mesh nano-metamaterial. 



complex nanoscale self-assembled metamaterial, such as 
the metallic gyroids [ll[ . The THM is composed of three 
wire helices oriented along the three orthogonal axes (as 
shown in the inset of figure^c)). Hence, as for the wire- 
mesh metamaterial, the THM shows a Drude behaviour 
for the effective electric permittivity: 



Xee 



iTuj 



u) + iTuj 



(11) 



where uj p = c yJp/{nR 2 a 2 (l - nR 2 /a 2 - L P /(2ttR))), V 
is a loss parameter and L is the self-inductance of the he- 
lical wires (see supplemental material and (lo| for more 
details). Unlike the wire-mesh, the effective magnetic 
response (xhh) of the THM is resonant and not al- 
ways equal to one. However, the most important char- 
acteristic of the THM is that it exhibits chiral com- 
ponents (keh and Khe), which in this case keh = 
—khe- Due to chirality, left-handed and right-handed 
circularly polarized waves propagating in a THM will 
experience different refractive indices, given by n± = 



1/ (yxS 



Xee^Xhh^) ± iK E x H {u)/co 



The electromagnetic response of a nanoscaled THM 
(i.e. the metallic components of THM are completely 
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penetrated by the incident electromagnetic fields) is de- 
rived by replacing all u terms in (|TT|) , xhh and keh 
with © (sec supplementary material) . Again, only the 24 
solution gives physical results. In figurc[5][a)-(b), we plot 
the analytical predictions for n_ and n + of a silver THM 
(with dimensions r w = 3.48nm, p = 20nm,R = 6.52nm 
and a = 40nm), along with numerical calculations that 
show remarkable agreement. The numerical calculation 
values were retrieved from S-parameter calculations, but 
using a different method as described in [2l|, due to the 
presence of chirality. As before, we change the geomet- 
rical parameters of the THM by varying the wire's ra- 
dius (r w ) and the host medium (e/j). The wavelength 
where Xee = for different structures is plotted in fig- 
ureOHc), where again the analytical (lines) and numerical 
(points) calculations show notable agreement. Remark- 
ably, as for the wire-mesh metamaterial, we observe a 
dramatic wavelength shift between the optical and mi- 
crowave models for extremely small geometrical changes 
in the structure. 

In summary, we have derived an analytical theory 
for nanoplasmonic metamaterials that accounts for the 
field penetration and plasmonic excitation of nano-scaled 
metallic geometries. This method opens a door to trans- 
fer established microwave analytical models to optical 
wavelengths. In fact, our theory accounts for metallic 
losses, which dominate the behaviour of optical metama- 
terials, any metallic dispersive behaviour and even in- 
ter sub-band transitions through the use of experimental 
data for the permittivity of metals. This method is a 
valuable tool to analytically predict the optical behaviour 
of any metamaterial and coupled optical nano-antennas, 
but most importantly complex self-assembled nanoplas- 
monic metamaterials, such as gyroidal structures. 
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